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Abstract. We investigate the influence of temperature and dissipation on the Landau-Zener 
transition probability in circuit QED. Dissipation is modelled by coupling the transmission 
line to a bath of harmonic oscillators, and the reduced density operator is treated within Bloch- 
Redfield theory. A phase-space representation allows an efficient numerical implementation 
of the resulting master equation. It provides reliable results which are valid even for rather low 
temperatures. We find that the spin-flip probability as a function of temperature and dissipation 
strength exhibits a non-monotonic behaviour. Our numerical results are complemented by 
Ch ' analytical solutions for zero temperature and for vanishing dissipation strength. 
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The demonstration of coherent quantum dynamics in superconducting flux and charge qubits 



OO , [1-3] represents a major step towards a solid-state implementation of a quantum computer. 

Manipulation and readout of a qubit can be achieved by a controlled interaction with an 



electromagnetic circuit. For charge qubits implemented with Cooper-pair boxes, this can 
be established either by a capacitive coupling to a transmission line [4, 5] or an oscillating 
circuit [6, 7] which can be modelled as a harmonic oscillator. The corresponding qubit- 
oscillator model also plays a role for the description of a flux qubit that couples inductively 
to an embracing dc SQUID [8]. These setups represent solid-state realizations of a two- 
level atom in an optical cavity [9, 10]. As compared to optical realizations, most solid-state 
implementations are characterised by a much larger ratio between qubit-oscillator coupling 
and oscillator linewidth [5]. 

Circuit QED experiments have already demonstrated quantum coherent dynamics [11], 
measurements with low backaction [7, 12], and the creation of entanglement between two 
qubits in a cavity [13, 14]. A further crucial prerequisite for a working quantum computer is 
quantum state preparation, i.e. the initialisation of the qubits. It has been suggested to achieve 
this goal by switching a control parameter through an avoided crossing with an intermediate 
velocity, such that the resulting Landau-Zener transition is neither in the adiabatic nor in the 
diabatic limit [15,16]. Then the avoided crossing effectively acts like a beam splitter. 
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The physics of Landau-Zener tunnelling is also of relevance for adiabatic quantum 
computation which relies on the time evolution of the ground state of a quantum system 
of a slowly time-dependent Hamiltonian [17]. In this scheme, a relevant source of errors 
are Landau-Zener transitions at avoided crossings between adiabatic energy levels. For an 
isolated two-level system, the corresponding transition probability has been derived in the 
classic works by Landau, Zener, and Stueckelberg [18-20]. When the qubit is coupled to 
a heat bath, this probability may change significantly. However, there exist also system- 
bath interactions for which the transition probability at zero temperature surprisingly is not 
affected at all by the coupling to a harmonic-oscillator bath [21] or a spin bath [22, 23]. For 
a heat bath with Ohmic spectral density, the transition probability at high temperatures is 
bath-independent as well [24-28]. The same holds true for the coupling to a classical noise 
source [29-31]. 

With this work we study finite temperature Landau-Zener transitions of a qubit that 
couples via a harmonic oscillator to a heat bath. For the qubit itself, this represents a case 
of a structured heat bath, because owing to its linearity, the oscillator plus the bath can be 
considered as an effective bath with a peaked spectral density [32-40] for which we provide 
results for the dissipative Landau-Zener problem for finite temperatures. In doing so, we 
restrain from eliminating this extra oscillator degree of freedom because in the present context, 
the oscillator dynamics is of experimental interest as well. The organisation of the work is as 
follows. In section [2l we introduce the qubit-oscillator-bath model and its treatment within 
Bloch-Redfield theory. Section [3] is devoted to Landau-Zener transitions for a thermal initial 
state, while in sectionlH we investigate dissipative transitions. Details of the derivation of the 
quantum master equation and its numerical treatment are deferred to the appendices. 

2. Model and master equation 

2.1. Landau-Zener dynamics in circuit QED 

Circuit QED involves a Cooper pair box that couples capacitively to a transmission line which 
is described as a harmonic oscillator [5]. The Cooper pair box is formed by a dc SQUID such 
that the effective Josephson energy E^ = £'j'cos(;r$/4>o) can be tuned via an external flux $, 
where $o denotes the flux quantum. We assume that 4> can be switched such that £; = hvt, 
V > 0, in a sufficiently large time interval. The capacitive energy ^Ec{N — Ng)'^ is determined 
by the number A^ of Cooper pairs on the island and the scaled gate voltage Ng. In the charging 
limit Ec ^ E], only the two states |A^ = 0) and |A'' = 1) determine the physics, and one defines 
the qubits states |T) °^ |0) + |1) ^^^ li) "^ |0) — |1)- Then, at the charge degeneracy point 
Ng = 1/2, the Hamiltonian reads 

J^s = -fi^<yz + Ji8(yx{a'^ +a)+naa'^a, (1) 

where the first term describes the qubit in pseudo spin notation with C7z|t,i) = ilT^i)- 
The second term refers to the coupling of the qubit to the fundamental mode of the 
transmission line which is described as a harmonic oscillator with the usual bosonic creation 
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Figure 1. Adiabatic energy levels of the qubit-oscillator Hamiltonian ([T]) as a function of the 
Josephson energy which is swept at constant velocity such that E] = hvt. The arrows mark the 
values of iij where the anti crossings are located. 

and annihilation operators a^ and a and the energy eigenstates |n) , n = 0, 1 , . . . , oo. Note that 
below we also consider the case of strong qubit-oscillator coupling for which a rotating- wave 
approximation for the coupling Hamiltonian is not justified. 

When the effective Josephson energy is switched from a large negative value to a large 
positive value, the system exhibits interesting quantum dynamics which can be qualitatively 
understood by computing the adiabatic energies as a function of time, see figure [U For 
low temperatures, we can assume that both the qubit and the oscillator are initially in their 
instantaneous ground state IT^O)- As time evolves, the system will adiabatically follow the 
state 1 1,0) until at time t = Q./v, an avoided crossing is reached. Then the system will evolve 
into the superposition a(v)|t,0) +/3(v)||, 1) with velocity-dependent probability amplitudes. 
This means that by adjusting the sweep velocity v, one can generate a single-photon state 
or qubit-oscillator entanglement [15]. For the time-evolution from t = —°° to f = oo, the 
corresponding bit-flip probability can be evaluated exactly and reads [15] 

-27tg^/v 



T-i 



1 



(2) 



Note that this generalisation of the Landau-Zener formula is also valid for large qubit- 
oscillator coupling ^ ^ ^ for which more than two levels are relevant and, thus, the scenario 
sketched above becomes more involved. 



2.2. Dissipative dynamics 

Dissipative effects in an electromagnetic circuit are characterised by an impedance Z(tt)) 
which, within a quantum mechanical description, can be modelled by coupling the circuit bi- 
linearly to its electromagnetic environment [41]. This provides the system-bath Hamiltonian 

[42-45] 

^ot = ^ + (a'^+a) Y,r^Ckibl + bk)+Y,^0)kblbk, (3) 

k k 
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where Mi = Mi + (a^ + aYY^fic^/lcOk is the Hamikonian of the qubit and the transmission 
line augmented by a counterterm, such that a frequency renormalisation due to the coupling 
to the bath is cancelled [42,43,45,46]. The second and third term describe the capacitive 
coupling of the transmission line to a bath of harmonic oscillators. The bath is fully 
characterised by its spectral density 

7((D)=£c|5((0fe-(0) = — W- wReZ((o). (4) 

k 

The second equality relates the system-bath model to classical circuit theory, where L and C 

are the specific inductance and capacitance of an effective transmission line that forms the 

dissipative environment. This relation can be established by comparing the resulting quantum 

Langevin equation with Kirchhoff's laws [41,47]. 

Since we are only interested in the dynamics of the qubit and the oscillator, all relevant 

information is contained in the reduced density operator p = trsPtotai which is obtained 

by tracing out the bath degrees of freedom. For weak system-bath coupling, the bath 

can be eliminated within Bloch-Redfield theory [48, 49] in the following way: Assuming 

that initially the bath is at thermal equilibrium and not correlated with the system, ptot °= 

p ^exp{ — Y,k^o^kbj^bk/k'BT), one derives within perturbation theory the quantum master 

equation 

^p = -^[^,p]-^°°dT{^(T)[e,[e(-T),p]]-fi^(T)[e,[e(-T),p]+]}, (s) 

where [A, 5]+ = AB + BA and [A, 5] = AB — BA denote the anti-commutator and the 
commutator, respectively, while the scaled "position" Q = a + a^ is the system operator to 
which the bath couples. The first term describes the unitary evolution generated by the system 
Hamikonian, while the second one captures the dissipative influence of the environment. 
The dissipative terms depend on the system dynamics through the Heisenberg operator 
Q{t) = e''^s?/7ige-i.^v7i jj^g i^^j]^ enters via the symmetric and the antisymmetric correlation 
functions 



1 f°° /' h03 \ 

^(0 = 2(['^(0.'^]+)eq= y^ A(OJ{(0)cOih[j^^yo^{m) 



(6) 



= -([<^(f),<^])eq = j d£0 7(£0)sin(£00 (7) 

of the effective bath operator t, =Y.kCk {b[ + bk). The Markovian master equation implies that 
the bath stays always close to equilibrium and that no system-bath correlations build up. 

For an explicit form of the master equation, we still need to evaluate the Heisenberg time 
evolution of the system operator Q{t). For a time-independent system Hamikonian, this can 
be done exactly by a decomposition into the energy eigenbasis. Then one obtains the standard 
Bloch-Redfield approach [48,49]. For periodically driven systems, the coherent dynamics is 
solved by the Floquet states which provide an appropriate basis [50]. The Hamikonian Qi, 
however, possesses a more general time-dependence and, thus, we have to resort to further 
approximations. If E] ~ hO., one could employ a rotating-wave approximation (RWA) for 
the cavity-qubit coupling, Ox{a + a'^) ^ o+a + O-a^ with o± — OciicTy. Within RWA, 
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the Hamiltonian ([T]) turns into the Jaynes-Cummings model for which the eigenvalues and 
eigenvectors are known [51]. Far off resonance, by contrast, i.e. for hO. ^ E] or fiQ. <^ Ej, 
an adiabatic approximation for either the qubit or the harmonic oscillator [52, 53] is helpful. 
For the present case of a Landau-Zener sweep, however, the Josephson energy assume all 
values from — oo to +°°, such that generally none of these approximations is well suited. 
Therefore, we resort to a weak-coupling approximation in the qubit-oscillator interaction 
g. The corresponding solution of the Heisenberg equations for the dimensionless position 



operator Q{t) is derived in the Appendix A and reads 

Q{t)=a'^c'^' + ac-'^'+g[k{t)o,-h{t)ay], (8) 

with the time-dependent functions Ic{t) = Q.[cos{Q.t) — cos{(Ojt)] / (Q?' — cof) and /s(f) = 
[coj sm{Q.t) — Q.cos{cO]t)]/{Q} — cof), where COj = Ej/h. In the derivation of this expression, 
we have neglected all terms of order (gt)^. Note that this approximation is used only for 
the evaluation of the dissipative kernels of the master equation ([5]), while for the solution 
of the master equation, the qubit-oscillator coupling is treated exactly. This means that we 
neglect dissipative terms of the order yg^ /Q? only, which is justified for a weak oscillator- 
bath coupling. 

For the further evaluation, we still need to specify the spectral density J{(o). Assuming 
that the environment is strictly ohmic and recalling that Q. = 1/v^ZC, we obtain from Eq. (Hj) 

•'('»' = 25a'" <" 

with the effective damping rate y. Inserting expressions ©-(E]) into Eq. ([5]), we arrive at the 



explicit quantum master equation (see again Appendix A ) 



^p = -^[^,p]-i^[e,[e,p]+]-|D„,[e,[e,p]] 

-4^^.p[e,[e,p]]-|i^cT[e,[c7„p]], do) 

with the operator Q = i//i[^, Q] = iQ.{a^ — a) and diffusion constants 

f ha \ 

'"'^2{co^+a^)^h^{v, + co,){v2,+a2y ^ ^ 

which refer to momentum diffusion and to cross-diffusion, respectively, while v„ = Ink^T /fi 
denotes the Matsubara frequencies and (Oc is a high frequency Drude cutoff. The prefactor of 
the last term contains an effective force 

which acts on the oscillator and depends on the state of the qubit. This means that the qubit 
dynamics influences the dissipative terms despite the fact that it couples to the bath only 
indirectly via the oscillator. This influence vanishes in the high-temperature limit, where in 
addition the diffusion coefficients become Dpp = Ik^T /hO. and Dj^p = 0, such that dissipative 
terms in Eq. ^ reduces to those of the well-known form derived in Ref. [54]. 
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The cross-diffusion diffusion term o^ D^p can be rather cumbersome due to its explicit 
dependence on the cutoff COc, which yields an ultraviolet divergence. Still it is possible to 
avoid its explicit evaluation and, thus, to render the cutoff obsolete with renormalisation 
arguments: The divergence and its regularization can be related to the ultraviolet divergence 
for the equilibrium momentum variance (2^)eq for which we find the relation 

^^xp=(e')eq-^'(e')eq, (14) 

where ( ■ ■ ■ )eq denotes the thermal average. It turns out that in the solution of the quantum 
master equation, D^p appears only in the combination (2^)eq = yDxp/^+Dpp (see Eq. (IB. 21) ). 
neglecting D^p is consistent with a weak-coupling approximation. Since for the harmonic 
oscillator the exact relation {Q^)eq = Dpp/Q} holds [55], the replacement D^p -^ provides 
the correct equilibrium expectation values, even in the limit of low temperatures. 

2.3. Solving the master equation in phase space 

The numerical solution of the quantum master equation (flOl ) requires an appropriate basis 
expansion. For the qubit, we choose the eigenstates of c^, i.e. 1 1) and 1 1) . The resulting matrix 
elements pij with iJ = ]A are operators in the Hilbert space of the harmonic oscillator. For 
these, at first sight, the natural basis is provided by the Fock states \n). With increasing 
temperature, however, good convergence is obtained only for a relatively large number 
of states. Then it is advantageous to transform the operators pij to a phase-space quasi 
distribution like the Wigner representation [56] which can be defined as 

The actual transformation can be accomplished by using Bopp operators [56]. Introducing for 
the qubit part the matrix notation 



M'j7(-»:.p) = ^/ dye'"'p,y(x--y,jr+-y). (15) 



'^^''''^^\ZZ^' *''' 



one obtains for the master equation ([101) the Wigner representation 

<9,W = ifhoW + ig[a„W]-f^<9p[a,,W]+-Fi^x[a,,W]-irF^^p[o:„W].(17) 

The Fokker-Planck-like operator 

ifho = ~^{pd, + xdp) + ydpp + yDppdl + jD^pd^p (18) 

governs the dissipative time-evolution of the harmonic oscillator, while the next three terms 
describe the coherent time evolution of the qubit and its coupling to the oscillator. The last 
term refers to the modification of the dissipative terms that stems from the qubit-oscillator 
coupling. 

The main advantage of this representation comes from the fact that the oscillator part ^o 
is formally identical with the Klein-Kramers equation [43, 57] for the classical dissipative 
oscillator, which allows one to adapt techniques for solving Fokker-Planck equations to 
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quantum master equations [50, 58, 59]. In particular, we will use the eigenfunctions 0„,/ (x, p) 
of =Sfho which obey the eigenvalue equation [60] 

^ho^nn'{x,p) = {nX+n*?i*)^„n>{x,p), n,n =0,1,2,..., (19) 

where A = —7/2 + 1(^^ — 7^/4) ^'^; see Appendix B. 1 The "ground state" ^00 is the Wigner 



representation of the density operator of the harmonic oscillator at thermal equilibrium. Thus, 
if the oscillator stays close to equilibrium, the decomposition of the density operator can 
be performed with only a few basis states — irrespective of the temperature. The resulting 



equations of motion for the coefficients can be found in Appendix B 



3. Landau-Zener tunnelling at finite temperature 

Thermal effects can modify the Landau-Zener transition probability ^ even in the absence 
of dissipation, i.e. for 7 = 0. Then the natural initial state is no longer the (initial) ground state 
1 1, 0) , but rather the canonical ensemble 

p(, = _oo) = |t)(t|® ^—, (20) 

where J^o = fid^a is the cavity Hamiltonian and ^ = ix[Qxg{—J^o/k^T)] the partition 
function at temperature T . Note that initially the qubit splitting is infinite, so that even at 
finite temperature, only the qubit ground state is populated. 

Since 7 = 0, the quantum master equation (flOl) reduces to the unitary von Neumann 
equation for the Hamiltonian ([I]), which is independent of the approximations made in the 
derivation of the master equation. Then the transition probability P^^i is the thermal average 
of the transition probabilities for the initial states 1 1 , «) , i.e. 

^T^T = Y^Pn P],n->] = Y^Pn P],n-.],m, (21) 

n n.m 

where p„ = Q-"^^/kYiT j a^ ^nd /'|,„^|,m = |(t,n|t/(oo;0)|t,m)p. It is worth noting that only 
terms with n > m contribute to the sum, due to the "no-go-up" theorem [22] which states 
that P'\^n->'\,m = for m> n. For the computation of the remaining probabilities P]_n^]_m we 
need to resort to an approximation: If all avoided crossings in the adiabatic qubit-oscillator 
spectrum (Fig. [T]) are well separated, one can treat the transitions at the avoided crossing as 
being independent of each other and compute the transition probabilities P]ji^]^m as joint 
probabilities [23]. In the vicinity of an avoided crossing, the qubit-oscillator Hamiltonian 
is described by the two-level system //tls = ^yto^ + k^^x with sweep velocity v and some 
level splitting A. The corresponding probability for a non-adiabatic transition is given by the 
standard Landau-Zener expression [18-20] 

w(A) = Q-''^^l^\ (22) 

For the Hamiltonian ([I]), the avoided crossings are formed only between states \],n) and 
I i, n + 1 ) at times t = ^Q./v with the level splitting 



A„ = IgvG^+l. (23) 
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Figure 2. Transition paths from the initial states |t,n), n = 0, 1,2, to the qubit state |t), i.e. 
those contributing to the probability P]^] in the individual-crossing approximation, which is 
vahd for weak qubit-oscillator coupling g<^Q.. 



Then the only paths that connect two qubit states |t) are sketched in Fig. |2l For the 
initial states |t,0) and |t, 1), the probabilities to end up in the qubit state |t) then become 



T,o^T 



w{2g) and Pj^ 



i^T 



w 



{2g)w{2g\/2), respectively. For oscillator states with n> \, 



two final oscillator states are possible. Assuming that interference terms do not play any role, 
we find the transition probability 

P^n-^] = w{2g^)w{2gV^iTl) + [1 -w{2gV^)][l ~w{2gVi[^)], (24) 

which formally also holds for n = 0, 1. Inserting (|24|) into (|2T1) . we obtain the transition 
probability P^^i = 1 —JlnPnP-\.n^-\^ where the sum can be identified as a geometric series. 
Evaluating this series, yields a central result, namely the Landau-Zener transition probability 
for finite temperature and weak qubit-oscillator coupling: 

n-i = fis'M - /(2^Vv) e-2-^Vv + [j(^2/^) _ /(2^2/v)] c^-s'/v ^ (25) 
where the temperature dependence is captured by the function 

m 



which for X > vanishes in the high-temperature limit, while f{x) = 1 for zero temperature. In 
the latter limit, expression (l25l) becomes identical to equation (|2l). The independent-crossing 
approximation is valid whenever the time between the anti-crossings, t = 2v/Q., exceeds 
the duration of an individual Landau-Zener transition, Tlz ^ 



/I7vmax{l,v^A27^} [61]. 
Inserting the explicit expression (|23l) for A„, this condition becomes Q. > ng. Thus, analytical 
result (|25l) hold only as long as oscillator states \n) with n > Q./g are not thermally occupied, 
i.e. for k^T < hQ?/g. Fortunately, in the range of current experimental interest, k^T < hQ. 
and ^ < ^, this condition is fulfilled. The numerical results shown in figures |3];a,b) confirm 
the results of the individual crossing approximation very well. Even for higher temperatures, 
k^T > hQ., we find very good agreement, provided that the qubit-oscillator coupling g is 
sufficiently weak, see figure[3tb). Moreover, the time evolution of the population (see figured 
below) confirms that the dynamics indeed discerns into two individual transitions. 

The temperature dependence of the spin-flip probability P^^i shown in figure[3]possesses 
an intriguing non-monotonic behaviour: For low temperatures, k^T < 0.2fiQ., the probability 
is almost temperature independent. With an increasing temperature, P]^[ first becomes larger, 
while it eventually converges to zero in the high-temperature limit. The temperature for which 
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Figure 3. Spin-flip probability P]^i as a function of the temperature for various sweep 
velocities and g — 0.04i2 (a) and for various qubit-oscillator coupling strengths and v — 0.01i2^ 
(b). The symbols mark the numerical solution of the Liouville-von Neumann equation 
(Eq. ( [TtI i for y — 0), while the solid lines refer to the individual-crossing result dZSl l. The 
dashed lines are a guide to the eye. 



2.5 



the the spin-flip probability assumes a maximum is essentially independent of the qubit- 
oscillator coupling g and increases with the sweep velocity. 

This behaviour can be understood in the individual-crossing scenario sketched in figure[2l 
For very low temperatures, initially only the state |t,0) is significantly populated and, thus, 
only the path sketched in figure l2t a) is relevant. For a temperature T ^ hQ./k^, also the initial 
state It, 1) becomes relevant. For this initial state, reaching the final state |t) requires two 
non-adiabatic transitions (see figure [2tb)) which enhances the probability to end up in state 
II). With further increasing higher temperature, states ||,n) with n>2 start to play a role. 
In the individual-crossing picture, each of these states can evolve into four possible finale 
states, two of which with spin up and two with spin down (see figure [2tc)). Consequently, the 
probability of reaching the spin state ||) is no longer enhanced. The relevance of oscillator 
states with n > 2 also qualitatively explains the fact that P^^i vanishes for high temperatures: 
In the adiabatic limit v ^ g^, the state | ],n) evolves via the state ||, n — 1) into the final state 
1 1, n — 2) . In the opposite limit of fast sweeping, the systems essentially remains in its initial 
state |n, t)- In both cases, the qubit will predominantly end up in state |t), which complies 
with our numerical results. 

The non-monotonic temperature dependence is in contrast to the behaviour found for 
Landau-Zener transitions of a qubit that is coupled to further spins, for which a monotonic 
temperature dependence has been conjectured [23]. The physical reason for this difference is 
that the spin coupled to the qubit possesses only one excited state. Then the paths sketched in 
figures [2b,c), which are responsible for the non-monotonic temperature dependence, do not 
exist. 
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4. Dissipative Landau-Zener transitions 

In the previous section we have studied the consequences of thermal excitations of the initial 
state for the transition probability ^ in the absence of an oscillator-bath coupling, i.e. for 
dissipation strength y =0. We next address the question how dissipation and decoherence 
modify Landau-Zener tunnelling. 

4.1. The zero-temperature limit 

For a heat bath at zero temperature, the exact solution of the dissipative Landau-Zener 
problem has been derived in recent works [21, 22]. Moreover, this limit generally is rather 
challenging for a master equation description of quantum dissipation [43,46]. Therefore the 
zero-temperature limit represents a natural test bed for our Bloch-Redfield master equation. 

Let us start with a brief summary of the derivation of the exact expression for the spin- 
flip probability /'j^j, as given in Ref. [22]. The central idea is to consider the cavity plus 
the bath as an effective bath that consists of "oo -|- 1" oscillators [32, 33, 35, 37^0]. Then the 
qubit-oscillator coupling operator Ox{a -\-a^) is replaced by a qubit-bath coupling of the type 
OxE/t'(^I' +^^')' where (a^ +^/t') denotes the normal coordinates of the effective bath. Thus 
the total Hamiltonian ([3]) can be written in terms of the spin-boson Hamiltonian [55] 

F (t\ 

^ot = ^——Oz-\-hgo^Y,^k'{dl,+dk')+Y,^^k'al,dk!. (27) 

^ V Id 

Moreover, the transformation to the normal coordinates provides the effective spectral density 
[32,62] 

/eff(a)) = g £c,,5(a) - CO,) = (^2_^2)2 + (y^)2 (28) 

with the effective dissipation strength 

a^-i^T (29) 

For the time-dependent Josephson energy E]{t) = hvt, this model defines a dissipative Landau- 
Zener problem for which at T = the spin-flip probability reads [22] 

P^^^ = l-Q-^''^^^k''l'l\ (30) 

Note that this formula is exact for any values of the coefficients g and Cj^i , provided that the 
systems starts dXt — — oo in its ground state. The sum Y^k'^lt can be expressed in terms of the 
spectral density (|28l) . such it becomes 

If. , . 1 



1 /■ 1 

E4 = — / dw7eff(w) = - 

k' s J ^ 

-2 




arctan — -f= 



(31) 



In the limit 7 -^ 0, we obtain Y,k' ^i ~ ^^ ^uch that the transition probability (1301) equals 
expression (O which is valid in the absence of the oscillator-bath coupling. For a small but 
finite dissipation strength y, Y^fj ^\t < 1 and, thus, the spin-flip probability becomes smaller 
when the oscillator is damped. 
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Figure 4. Comparison of the spin-flip probability P]^i at T = obtained with Bloch-Redfield 
theory (symbols) and the exact result, Eqs. dSOl ) and (ISTT l. as a function of the dissipation 
strength for qubit-oscillator coupling strength g = 0.04i2 and various sweep velocities. The 
probability has been normalised to the corresponding value in the absence of dissipation, 
P^^l (7=0). The deviation from the exact result is always below 1 %. 



We now use this exact result as a test for the master equation (iTOl) . It is worth emphasising 
that both problems are not fully equivalent, because the preparations are different: At zero 
temperature, the initial condition (|20l) for the master equation reads p{t = —°°) = 1 1, 0) (t, 0| , 
i.e. the composed qubit-oscillator-bath system starts in the pure state ||,0) ® |0,0, . . .), where 
the latter factor refers to the bath. In the analytical treatment sketched in the preceding 
paragraph, by contrast, the initial condition is the ground state of the total Hamiltonian (l27l) . 
1 1) ® |0, 0, 0, . . .). Since a finite oscillator-bath coupling induces system-bath correlations [45], 
the oscillator is generally entangled with the bath. Nevertheless, in the present context, the 
difference should be minor, because we consider a time evolution that starts att = —0°, such 
that the oscillator-bath setting can evolve into its ground state before Landau-Zener tunnelling 
occurs. 

In figure HI we compare the transition probabilities obtained with the quantum master 
equation (fTOl) with the corresponding exact analytical result (l30l) . For the system parameters 
used below, we find that Bloch-Redfield theory predicts even at zero temperature the exact 
results with an error of less than one percent. Since the quality of a Markovian quantum 
master equation generally improves with increasing temperature, the results presented below 
are rather reliable. 



4.2. Thermal excitations and dissipative transitions 

We next turn to the generic situation in which both thermal excitations of the initial state 
and dissipative transitions play a role, i.e. we consider the situation of finite temperatures 
and finite dissipation strength. The resulting spin-flip probabilities P]^i for three different 
sweep velocities are shown in figure [51 For small temperatures, k^T < 0.2^Q, we find 
that dissipation slightly reduces the spin-flip probability P]^[. This is consistent with the 
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Figure 5. Damping effects on the asymptotic transition probability P\^i as a function of the 
temperature for qubit-oscillator coupling g — 0.04i2, various sweep velocities v and damping 
strengths 7. The symbols depict numerical results obtained with the quantum master equation 
( fTTl ). while the solid lines mark the results from the independent-crossing approximation ( 
for 7=0. 



behaviour at zero temperature, discussed in section 14.11 Once the thermal energy is of the 
order liQ, the opposite is true: dissipation supports transitions to the final ground state ||) 
and, thus, P\^i increases. This tendency is most pronounced for the intermediate sweep 
velocity chosen in figure Ob). In this regime we again find a non-monotonic temperature 
dependence of the transition probabilities. 

In order to reveal the role of dissipative decays, we focus on the population dynamics 
for an intermediate temperature k^T ^ O.ShQ., where the transition probabilities are already 
significantly influenced by thermal excitations; cf. figure [5l Nevertheless, the temperature is 
still sufficiently low, such that only the oscillator states |0) and 1 1) are relevant. Figure[6] shows 
the time evolution of the population of the states |t,0) and |t, !)• Obviously, the populations 
change considerably at the avoided crossings of the qubit-oscillator spectrum, so that the 
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Figure 6. Population dynamics of the states |t,0) (a) and |t, 1) (b) during the Landau-Zener 
sweep. Dissipation strength and temperature are 7 = 10^^i2 and k^T — O.ShQ., respectively. 
The gray lines are the corresponding results in the zero-friction limit 7=0, for which the 
probabilities converge to the values indicated by the dashed lines. The arrows mark the 
positions of the avoided crossings in the qubit oscillator spectrum; cf. figure[T] 



dynamics discerns into three stages. 

For t < —0./v, the system remains in the canonical state. When at time t = —Q./v the 
anti-crossing between the states |t,l) and ||,0) is reached, the population of the former 
state drops due to an adiabatic transition to the latter. As a consequence, the oscillator is 
no longer at thermal equilibrium and, consequently, we observe thermal excitations from 
1 1,0) to It,!)- When at time t = Q./v, the second set of anti-crossings is reached (see 
figure [2ta,b)), both states undergo an individual Landau-Zener transition after which the 
populations converge in an oscillatory manner towards their final value. At the final stage, 
the oscillator populations thermalize, while those of the qubit remains practically unchanged. 
The latter sounds counter-intuitive because the physical system is dissipative. Nevertheless 
this has a physical explanation: The qubit experiences an effective heat bath with a spectral 
density sharply peaked at the oscillator frequency Q.. Thus for large times, t ^ ^/v, the 
spectral density at the qubit splitting hvt vanishes and, consequently, the qubit is effectively 
decoupled from the bath. 

The population dynamics provides evidence that for weak dissipation and narrowly 
avoided crossings, the dynamics consists of individual Landau-Zener transitions, while 
dissipation takes place mainly in-between the transitions. This behaviour resembles the one 
occurring in nanomagnets [63,64]. 
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5. Conclusions 

We have investigated the influence of finite temperature, decoherence, and dissipation on 
Landau-Zener transitions of a two-level system (qubit) that is coupled via a harmonic 
oscillator to a heat bath. In particular, we have focussed on a recent solid-state realization 
of this model, namely the so-called circuit QED for which Landau-Zener sweeps can be 
performed by switching the effective Josephson energy of the Cooper pair box. The adiabatic 
spectrum of this system consists of a sequence of exact and avoided crossings, where 
for strong qubit-oscillator coupling, the latter may even overlap. Therefore, the resulting 
quantum dynamics is more complex than in the "standard" two-level Landau-Zener problem. 
Moreover, this qubit-oscillator-bath model is equivalent to coupling the qubit to a bath with 
peaked spectral density. 

Dissipation has been modelled by coupling the oscillator to an Ohmic environment which 
we integrated out within a Bloch-Redfield approach. We solved the resulting master equation 
numerically after a transformation to Wigner representation followed by a decomposition into 
proper basis functions. The comparison with results for the exactly solvable zero-temperature 
limit, demonstrated that our approach provides reliable results, even though this limit is known 
to be rather challenging for a Markovian master equation. 

For vanishing dissipation strength, the temperature enters only via initial thermal 
excitations of the oscillator. Most interestingly, we found for this case that the spin-flip 
probability exhibits a non-monotonic temperature dependence. For a sufficiently small qubit- 
oscillator coupling, this can be understood within the approximation of individual Landau- 
Zener crossings. This picture reveals the special role played by the first excited oscillator state. 
As compared to any other state, this state is more likely to induce a spin-flip. At intermediate 
temperatures, the first excited oscillator state possesses a relatively high influence, which 
leads to the observed non-monotonic behaviour. When dissipative decays become relevant 
as well, transitions to the final (adiabatic) ground state of the qubit become more likely. 
Nevertheless, in some small regions of parameter space, we find the opposite, namely that 
the probability to find the qubit in the excited state is slightly enhanced. This is at first sight 
counter-intuitive, but can be understood from the fact that the qubit effectively experiences 
a bath with a peaked spectral density. Therefore, dissipative qubit transitions can occur only 
during the short lapses of time in which the adiabatic energy splitting of the qubit is of the 
order of the oscillator frequency. With an increasing oscillator-bath coupling, the peak in 
the spectral density becomes broader and, thus, the more intuitive tendency towards the final 
ground state starts to dominate. 

Finally, our results provide evidence that the recently derived zero-temperature results 
hold true also for finite temperatures, provided that the thermal energy does not exceed a 
value of roughly 20 percent of the oscillator's energy quantum. This means that experimental 
quantum state preparation schemes that rely on Landau-Zener transitions in the zero- 
temperature limit, are feasible already when the oscillator is initially in its ground state, while 
the low-frequency modes of the bath may nevertheless be thermally excited. 
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Appendix A. Derivation of the quantum master equation 

In this Appendix we outline the derivation of the master equation (flOl) starting from the general 
Bloch-Redfield expression ([5]). 

Appendix A. 1. Heisenberg coupling operator 

The essential step is the solution of the Heisenberg equation of motion for the scaled position 
operator Q = a-\-a^ oi the oscillator, which will rely on approximations. In doing so, we even 
address a slightly more general qubit Hamiltonian outside the charge degeneracy point, i.e. 
we also consider the charging energy jEqOx, such that the Rabi Hamiltonian ([T]) becomes 

J^s = -^Oz-^Ox + hQ.a^a + hgOjc{a^ + a). (A.l) 

For convenience, we write this Hamiltonian in the eigenbasis of the qubit: 

hA 
J^s = — —d^ - hQa' a + hg{cos do^ - smdOx){a' + a) , (A.2) 

where 



hA = \/E]^+Ec^, e = arctan — (A.3) 

V Ec 

denote the energy splitting and the coupling angle, respectively, of the qubit. The 

corresponding Heisenberg equations become 

Q = -Q}Q-ga{cosea.-smeax), (A.4) 

Oz= -2gsm0OyQ, (A.5) 

dx = {A- 2g cos 0)Oy, (A.6) 

dy= -{A-2gcosG)Ox + 2gsm0o.x, (A.7) 

which are non-linear due to the qubit-oscillator coupling and, thus, cannot be solved directly. 
We are only interested in the time evolution of the coordinate Q of the oscillator with 
eigenfrequency ^. Since in typical circuit QED experiments ^ ^ ^, we treat the coupling as 
a perturbation. In the absence of the coupling, the time evolution of the qubit operators reads 

Ox{t) = OxCos{At) +OySm{At), (A.8) 

a,(0 = o,. (A.9) 

Inserting this into the equation of motion (IA.4I) for the oscillator coordinate, it becomes 
evident that the qubit entails on the oscillator the time-dependent force 

F{t) = -gQ.lozCos0- [axCos{At) + OySm{At)] sinol. (A.IO) 
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To first order in g, the solution of equation (IA.4I ) reads 






Q{t)^ae''"+a'^ &-''"+ dt' G{t -t')F{t'), (A.ll) 

Jo 

where 

GW = ^9W (A.12) 

denotes the retarded Green function of the classical dissipative harmonic oscillator. 
Evaluating the integral we finally obtain the Heisenberg operator 

Q{t)=aQ'^'+a^Q-''^' -gh{0,a;t)o^co^e 

+ g[h{Q.,^■,t)o^ + h{Q.,^■J)Oy] sine, (A.IB) 

with the functions 

, , , aco^iat) — aco^ibt) 

IciaMt) = 2 ,2 , (^-14) 

a^ — b^ 



hiaMt) = , ,. • (A. 15) 



bsm{at) —asm{bt 
a^-b^ 
Expression (IA.13I ) allows one to explicitly evaluate the Bloch-Redfield equation ([5]). 



Appendix A.2. Ohmic spectral density 

In circuit QED, the environment of the qubit and the transmission line is formed by electric 
circuits and, thus, can be characterised by an effective impedance. In most cases, this 
impedance is dominated by an Ohmic resistor, which corresponds to the Ohmic spectral 
density ^ of the bath. Then the time integration in the Bloch-Redfield equation ([5]) can 
be evaluated. The antisymmetric correlation function ([7]) is then given by 

^ ^ dTft)sin(ft)T) = --^— 5(t), (A.16) 



27tQ.Jo ^ ' 2a dT 

such that the last term of equation ^ becomes {iY/4Q.)Q. The remaining time integrals are 
of the type 

J^ dT^(T)cOs(aT) = -COth(^^j (A.17) 

/ dT^(T)sm(aT) = ^= — -^ > -—^ — — ^ , (A.18) 

where we introduced the Matsubara frequencies v„ = Innk^T /h. In the second integral, an 
ultraviolet divergence has been regularised by a Drude cutoff e^®'®'= [44]. 

Appendix B. Basis expansion 

In this appendix we outline the diagonalisation of the oscillator Liouvillian in Wigner 
representation, =^ho5 whose eigenvectors are used as a basis set for the numerical treatment. 
Since the operator ^^o^ apart from the cross-diffusion D^p, is of the same form as the Fokker- 
Planck operator of the corresponding problem for classical Brownian motion, we can make 
use of an idea put forward by Titulaer [57, 60] and generalise it along the lines of Ref. [50]. 
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Appendix B.l. Diagonalisation of the oscillator Liouvillian 

By solving the characteristic functions of the partial differential equation ^ — =2ho0, one finds 
the operators 

r^ = d^ + -dp, (B.l) 






^,_^,{{Q'US,-^D„S,+.--p), (B.2) 

which commute with d /dt — =^ho and, thus, map any solution of the Liouville equation to 
a further solution. For notational convenience, we have introduced the eigenvalues of the 
classical equation of motion of the dissipative harmonic oscillator. 



A = -|+iya2-^, (B.3) 

and A*. It will turn out that the equilibrium expectation value of the dimensionless oscillator 
coordinate Q becomes (2^)eq = Dpp + yD^pjO.. Since moreover, the diffusion coefficient 
D^p appears in all results only in this combination, the introduction of (2^)eq turns out to be 
convenient as well; see also the discussion after equation (fT3l) . The operators (IB. II) and (IB. 21) 
fulfil the commutation relations 

[r_,r+] = l, [r_,r;] = [rl,r+]=0, [ifho,r±] = ±Ar± (B.4) 

and allow one to write the Liouvillian in the diagonal form 

^ho = A r+r_ + A V;rl, (B.5) 

where the symbol * denotes complex conjugationla Formally we have reduced the eigenvalue 
problem for the Liouvillian to that of two uncoupled harmonic oscillators, so that the 
eigenvalues obviously read 



A„„, = An + AV = -(n + n')| + i(n-n')Y^2_I_^ (B.6) 

where n, n' = 0, 1, 2, . . . The corresponding eigenfunctions 0„,/ can be constructed by applying 
the raising operators r^ and r\ on the "ground state". The latter is the stationary state defined 
by the relation ^ho^oo = whose solution is the Gaussian 

Since ^ho 7^ -^ho' ^^ eigenfunctions 0„„/ are not mutually orthogonal, so that we need 
to compute the left eigenvectors ^,„j/ of ^ho as well. Repeating the calculation from above for 
°^ho' ^^ ^'^^ ^^ ^^^^ ground state ^oo = 1, so that we obtain the eigenfunctions 

<^™' = ^(r+)"(r;)">oo (B.8) 

n\n \ 

^nn' = {r+nrir' 1 (B.9) 

t We do not consider the overdamped limit 7 > 2 A in which A becomes real and, thus, the notation needs to be 
modified. 
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which fulfil the ortho-normalisation relation 

/ docdp ^mm' <^nn' = ^nm ^n'm' ■ (B . 1 0) 

Appendix B.2. Expansion of the entire Liouvillian 

Besides the already diagonalized oscillator Liouvillian, the total Liouvillian (flTI ) for the qubit 
plus the oscillator contains also the operators 

X = -r^-r*_ + Kr+ + K*r% ^ = (Q^)^<i ^2 _q2 ^ (^-^l) 



r+-r* 



dp=-i — , + (B.12) 

^ 2^1-72/4^2 

Then, the basis decomposition of the Wigner function, 

W = ^C„y^„„'; C„y = I "!|" "j" 1, (B.13) 

nn' \ '^n.n' ^n,n' J 

obeys the equation of motion 

Cn,n' = {nX+n A*)C„y +i— [a2,C„y] -i^[crt-,C„+i,„/ + C„,„/+i] 

a/1 - 72/4^2 

g 

v/i - 72/402 ^ ' ' ^ 

- [c7x,«'C„y_i-nC„_iy], (B.14) 



'1 - 72/4^2 

which is a set of 4A^2 coupled linear ordinary differential equations. 

Appendix B.3. Computation of expectation values 

The expectation value of an operator can be performed directly in the basis of the 
eigenfunctions without back transformation to the operator representation of the density 
operator. For an observable 

B = S0O{x,p), (B.15) 

for which S and O refer to qubit and oscillator variables, respectively, the expectation value 

(5)=tr(5p) =Y,S,j fdxdpO{x,p)Y,cl„,^nn'ix,p) (B.16) 

ij nn' 

can be expressed in terms of the operators (|B.1I) and (IB. 21) such that 

(5)=££5,,O„,,c;/;„,(0, (B.17) 

ij nn' 

where 0{x,p) is the corresponding operator in phase-space [56] and 

Onn' = / djcd;? 0{X, p) (j)nn' (x, p) ■ (B. 1 8) 
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If one is only interested in the behaviour of the qubit, i.e. for 0=1, the expectation value 
becomes 

{B) = {S)=Y,SijC%{t). (B.19) 

ij 

This implies that any information about the qubit state is already contained in the four 
coefficients c^{t), which represents a particular advantage in the present decomposition. 
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